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The Equivalence of Semidefinite Relaxation 
MIMO Detectors for Higher-Order QAM 
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Abstract 

In multi-input-multi-output (MIMO) detection, semidefinite relaxation (SDR) has been shown to be 
an efficient high-performance approach. Developed initially for BPSK and QPSK, SDR has been found to 
be capable of providing near-optimal performance (for those constellations). This has stimulated a number 
of recent research endeavors that aim to apply SDR to the high-order QAM cases. These independently 
developed SDRs are different in concept and structure, and presently no serious analysis has been given 
to compare these methods. This paper analyzes the relationship of three such SDR methods, namely 
the polynomial-inspired SDR (PI-SDR) by Wiesel et al, the bound-constrained SDR (BC-SDR) by 
Sidiropoulos and Luo, and the virtually-antipodal SDR (VA-SDR) by Mao et al. The result that we have 
proven is somehow unexpected: the three SDRs are equivalent. Simply speaking, we show that solving 
any one SDR is equivalent to solving the other SDRs. This paper also discusses some implications arising 
from the SDR equivalence, and provides simulation results to verify our theoretical findings. 
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I. Introduction 

Multiple-input-multiple-output (MIMO) detection using semidefinite relaxation (SDR) [1]— [14] has 
recently received increasing attention. Being a symbol-constellation dependent technique, SDR has been 
shown to provide considerably better symbol error performance than some suboptimal MIMO detectors 
such as the linear and decision-feedback receivers. SDR is not an optimal approach from a maximum- 
likelihood (ML) perspective, but it guarantees a polynomial-time complexity in the problem dimension. 
In comparison, the currently best known methods for optimal ML MIMO detection, namely sphere 
decoding [15], [16] do not have such a guarantee [17]Q 

SDR was first proposed for the BPSK constellation [1], [2], and the same idea can easily be carried 
forward to the QPSK constellation (or 4-QAM) [6], [7]. For BPSK and QPSK, simulation results have 
indicated that SDR can provide near-optimal bit error performance. This intriguing finding has stimulated 
a number of works. Theoretically, it is shown recently [13] that SDR can achieve the full receive diversity 
in the BPSK case. An equally interesting but totally different analysis is given in [14], where the SDR 
approximation gap is examined using random matrix theory. Apart from theoretical analysis, there has 
been interest in various aspects such as fast implementations [19], [20] and soft-in-soft-out MIMO 
detection [7]. 

But what attracts more attention in SDR is possibly its extension to more general symbol constellations, 
especially the higher-order QAM. For higher-order PSK, readers are referred to [3]. For higher-order QAM 
which is the focus of this paper, the first endeavor is by Wiesel et al. [9], who proposed a polynomial- 
inspired SDR (PI-SDR) method for 16-QAM. In that work the authors also showed that PI-SDR is a 
bidual of the optimal ML MIMO detector (or achieves an optimal Lagrangian dual lower bound of the 
ML). The drawback of PI-SDR is that its extension to larger QAM sizes would be increasingly complex 
to handle. Later, Sidiropoulos and Luo proposed a bound-constrained SDR (BC-SDR) method [9] for 
any QAM. BC-SDR aims at simplicity and appears to be less sophisticated than PI-SDR. For instance, 
the BC-SDR problem structure is virtually the same for any QAM size. The simple, special structures of 
BC-SDR make fast implementations [21] particularly favorable. More recently, Mao et al. [10] developed 
a virtually-antipodal SDR (VA-SDR) method for any ^-QAM (where q > 1). As its name implies, VA- 
SDR has a strong connection to the SDR used in BPSK/QPSK. VA-SDR is structurally less complex 
than PI-SDR, but involves more optimization variables than BC-SDR. 

'As a short aside, the complexity limitation of optimal sphere decoding has recently motivated interest in some suboptimal 
but reduced-complexity variants; e.g., the Fano decoder [16] and the fixed-complexity sphere decoder [18]. 
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It is worthwhile to mention two more recent developments. Mobasher et al. formulated a class of SDR 
problems that is applicable to any kind of symbol constellations [12]. As a price for their generality, 
Mobasher's formulations are considerably more complex than the other SDRs. This translates into a 
higher computational requirement. Yang et al. [11] proposed a tightened version of BC-SDR for the 16- 
QAM case. Interestingly, they showed that their tightened BC-SDR can provide a better approximation 
than the 16-QAM PI-SDR. 

While a number of SDR methods have been proposed for higher-order QAM, their comparisons have 
not been seriously considered at present. This paper focuses on analyzing the relationship of the PI-SDR, 
BC-SDR, and VA-SDR methods. We obtain a result that is intuitively not so obvious: PI-SDR, BC-SDR, 
and VA-SDR are actually equivalent, despite the fact that they exhibit rather different problem structures. 
Specifically, our analysis reveals that 

1) For 16-QAM and 64-QAM, PI-SDR is equivalent to BC-SDR. 

2) For any 4 g -QAM, VA-SDR is equivalent to BC-SDR. 

The SDR equivalence result to be established is quite general, and it can be applied to problems other 
than MIMO detection. For example, the SDR equivalence is perfectly applicable to the problem of blind 
ML detection of orthogonal space-time block codes [22] where the objective function structure is different 
from that in MIMO detection. 

We attempt to give insights into the SDR equivalence result by organizing this paper in the following 
way. We use the relatively simple case of 16-QAM to provide the problem statement in Section [TTJ and 
to review the three SDR methods in Section |III] This is followed by Section |lVl where we describe 
the 16-QAM SDR equivalence, illustrate the ideas behind, and discuss its implications. Then, Section |Vl 
addresses the more complex case where the QAM size is further increased. Some simulation results are 
provided in Section [VT] for verifying our claim and for demonstrating the SDR performance. 

II. Problem Statement 

We consider a standard MIMO system model, namely that of spatial multiplexing (or V-BLAST) over 
a frequency-flat channel [23], [24]. Supposing that the system has N transmitter antennas and M receiver 
antennas, the received signal can be described by the following formula: 

y = Hs + i>. (l) 

Here, H G £,MxN j s ^ jyjijyjQ channel, y G C M is the received signal vector, u is a noise vector 
assumed to be zero-mean circular white Gaussian, and s G S N is the transmitted symbol vector where 
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S C C denotes the symbol constellation set. For the 16-QAM constellation, S is given by 

S = { s = s R + jsi | s R , si e {±1, ±3} }. 

It should be stressed that in many digital communication problems, the signal models can be formulated 
in the same form as in CO); for example, the CDMA multiuser problem. Readers can find more examples 
in the literature such as [15], [16]. Thus, detection techniques for (Q]), or simply MIMO detection, have 
wide applicability. 

It would be convenient to reformulate the complex-valued model in £lj to a real-valued one. Let 
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M = 2M, and N = 2N. Eq. © is equivalent to 

y = Hs + v (2) 

where s £ {±1, ±3}^ and v is defined in the same way as y. The ML detection problem for the MIMO 
model in $2} is shown to be an optimization 

min lly — Hsll 2 , (3) 

sg{±l,±3}« 

in which the globally optimal solution serves as the ML decision. Note that || • || in (f3]) stands for the vector 
2-norm. ML detection is known to provide superior detection performance, but the major challenge lies 
in solving ([3]) which is a computationally hard problem. Presently, the best known optimal solver for © 
is sphere decoding [15], [16]. While sphere decoding has been empirically found to be computationally 
very fast for small to moderate problem sizes (say, for N < 20 for 16-QAM), it has been revealed [17] 
that the sphere decoding complexity would be prohibitive for large N and/or low SNRs. 

III. Review of Three 16-QAM SDR Detectors 

SDR is a suboptimal approach to ML, using a class of polynomial-time solvable convex optimization 
problems known as semidefinite programs. In this section, we review three SDR methods for the 16-QAM 
constellation, namely PI-SDR [8], BC-SDR [9], and VA-SDR [10]. (Their extensions beyond 16-QAM 
will be considered later in the paper.) 
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A. Polynomial Inspired SDR 

PI-SDR was the first application of the SDR principle [25] to 16-QAM ML detection, to the best of 
our knowledge. To present its idea, consider a reformulation of the ML problem in ([3]) 

min tr(H T HS) - 2s T H T y + llyll 2 
SeS^seR" (4) 

s.t. S = ss T , Sue {1,9}, i = l,...,N 
where S is a slack variable, is the set of N x N real symmetric matrices, SV, denotes the (i, j)th 
element of S, and tr(-) is the trace operator. PI-SDR was inspired by the fact that 

u G {1, 9} (u - l)(u - 9) = u 2 - Wu + 9 = 0. 

By turning the constraints Su G {1, 9} to a polynomial form, Problem © is further reformulated as 

min tr(H T HS) - 2s T H T y + ||y|| 2 

S,s,U,u 

s.t. S = ss T , U = uu T (5) 
d(S) = u, d(U) - lOu + 91jv = 
where d : R NxN -> R N is the diagonal operator (i.e., d(S) = [ S n , S NN ] T ), and 1^ is the N- 
dimensional all-one vector. The reformulated ML problem in © is still hard, where the difficulty lies in 
the nonconvex constraints S = ss T and U = uu T which restrict S and U to be of rank 1. 
In PI-SDR, we relax the polynomial ML formulation in ([5]) to 

min tr(H T HS) - 2s T H T y + ||y|| 2 

s.t. Shss T , U^uu T (6) 
d(S) = u, d(U) - 10U + 91AT = 
where A ^ B means that A — B is positive semidefinite (PSD). The idea is to replace the hard constraint 
S = ss T by a convex constraint S y ss T , and similarly to U = uu T . There are two basic advantages 
with such a relaxation. First, Problem ([6]), or the PI-SDR problem is a semidefinite program (SDP) which 
is convex and does not suffer from local minima. Second, as an SDP the PI-SDR problem can be solved 
by available interior-point methods [26], [27] in a polynomial-time fashion. 

Once we solve the PI-SDR problem in (O, we can make a symbol decision by simple rounding of the 
PI-SDR solution associated with s. A better alternative to this simple rounding is the Gaussian randomized 
rounding; see [2], [8], [9] for the details. 

It has been shown that PI-SDR is a bidual of ML, a desirable property from a Lagrangian duality 
viewpoint. Specifically, © achieves the optimal Lagrangian dual lower bound of the polynomial ML 
formulation in d5). 
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B. Bound constrained SDR 

BC-SDR is possibly the simplest among the various 16-QAM SDR methods. It relaxes the ML problem 
in © to an SDP 

min tr(H r HS) - 2s T H T y + llyll 2 

V J 11*11 (7) 

s.t. S^ss T , l<S li <9, i = l,...,N, 
where the original constraint S = ss T is replaced by the PSD constraint S y ss T (as in PI-SDR), and 
the discrete set {1,9} is relaxed to an interval [1,9]. 

The BC-SDR problem in ([/J exhibits particularly simple SDP problem structure. This has enabled us 
to develop a specialized interior-point algorithm for |7J that runs many times faster than some general- 
purpose interior-point software [21]. 

While PI-SDR is bidual of ML, it is not known if BC-SDR has such a desirable property. 

C. Virtually Antipodal SDR 

VA-SDR was proposed by Mao et al. [loj^l who also suggested a MIMO detector that combines VA- 
SDR and multistage decision-aided cancellation. Here we are interested only in VA-SDR. The idea stems 
from the fact that 

se {±l,±3}^s = 6i + 26 2 , 6i,6 2 e{±l}. 
Hence, the 16-QAM ML problem can be re-expressed in a virtually antipodal form 

min ||y-H(bi + 2b 2 )|| 2 = min lly - HWbll 2 , (8) 

bi.b^il}* ' b€{±l} 2 " 

where we denote 

W=[I2I], b=[bj blT- 

By applying the same SDR as in BPSK/QPSK constellations, VA-SDR is obtained: 

min tr(W T H T HWB) - 2b T W T H T y + llyll 2 

I ) 11*11 (9) 

s.t. B y bb , B u = l, i = l,...,2N. 
In terms of problem structure, VA-SDR is exactly the same as the SDR for BPSK/QPSK. Hence, VA-SDR 
can be implemented by directly applying interior-point algorithms designed for BPSK/QPSK SDR [20], 
[26]. 



2 In fact, an earlier work by Steingrimsson et al, [7] was close to finding VA-SDR. In that paper, a symbol is considered as 
a linear transformation of some bit symbols, which is exactly how VA-SDR works. But we should emphasize that it was Mao 
et al. [10] who first described the use of VA-SDR for higher-order QAM and put the method to the test. 
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It can be shown that VA-SDR is also a bidual of ML. Specifically, (O achieves the optimal Lagrangian 
dual lower bound of the virtually antipodal ML formulation in ([8]); see the literature such as [28, 
Section 5.15 and Exercise 5.39] for the proof. We however should point out that the biduality of VA-SDR 
and PI-SDR does not translate into equivalence of the two SDRs. The reason is that the two SDRs are 
obtained from two different ML formulations, which may be viewed as different optimization problems 
that have the same optimal value. The two ML formulations could yield different dual optimal values 
(and consequently different SDR optimal values), unless the two formulations have strong duality which 
is unlikely in this problem. 

D. Comparisons by Simulations 

Let us use simulations to compare the above described SDR methods, before proceeding to analyzing 
them. The simulation setting follows that of a standard MIMO system, where the channel matrix H is 
i.i.d. complex circular Gaussian distributed with zero mean and unit variance. The MIMO system size is 
(M,N) = (8,8). For PI-SDR and BC-SDR, we employ the simple rounding procedure; i.e., if s* is the 
PI-SDR/BC-SDR solution associated with s, then 

s = dec(s*) 

is the detected symbol vector where dec(-) is the elementwise decision function for the discrete set 
{±1, ±3}. For VA-SDR, there are two possible ways of doing simple rounding. Let b* = [ (h\) T (b^) T ] T 
(b*, b£ G R^) be the VA-SDR solution associated with b. We can detect s either by 

s = sgn(bj) + 2sgn(b^) (10) 

where sgn(-) denotes the elementwise sign function, or by 

s = dec(b^ + 2b* 2 ) = dec(Wb*). (11) 

We call (fTOl and (fTTT > simple rounding I and simple rounding II, respectively. 

The simulated symbol error performance of the three SDRs is given in Fig. [TJ In the figure the SNR 
is defined as the received signal-to-noise ratio per QAM symbol; i.e., E{||hjSj|| 2 }/E{||i>|| 2 }. One can 
see that for VA-SDR, simple rounding II gives better performance than simple rounding I. But, more 
importantly, the performance of PI-SDR, BC-SDR, and VA-SDR (with simple rounding II) is identical. 
From this observation it is reasonable to suspect that there are strong connections between the three 
SDRs. 
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Fig. 1. Comparison of symbol error rates of PI-SDR, BC-SDR, and VA-SDR in an 8 x 8 16-QAM system. 



IV. Equivalence of the Three 16-QAM SDR Detectors 

In this section we prove the equivalence of PI-SDR, BC-SDR, and VA-SDR in the 16-QAM case. 
In the first subsection, the equivalence result and its implications will be described. Then, the analysis 
leading to the equivalence result will be shown in details in the second and third subsections. 

A. Main Result and Implications 

The three SDRs can be represented by a unified expression 

min /(S,s) (12) 

where 

/(S, s) = tr(H T HS) - 2s T H T y + ||y || 2 

is the objective function, and T is the feasible set, the definition of which depends on the SDR method 
employed. For BC-SDR, the feasible set is defined as 

^BC-SDR = { (S, s) | S h ss T , 1 N X d(S) H 91 N } (13) 

(We adopt the standard notation that and > y' mean elementwise inequalities, when applied on vectors). 
For PI-SDR, the feasible set is characterized as 

^PI-SDR = { (S, S) | (U, U, S, S) G W P |_SDR } (14) 

Wpi_sdr = { (U,u,S,s) I U y uu T ,S y ss T ,d(S) = u,d(U) - 10u + 91iv = }, (15) 
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and for VA-SDR, 

^VA-SDR = { (S, S) = (WBW T , Wb) I (B, b) G £vA-SDR I (16) 

fivA-SDR = { (B, b) | B y bb T , d(B) = 1 2N }. (17) 
Essentially, the equivalence of the three SDRs lies in the feasible set: 

Theorem 1 The feasible sets of the three 16-QAM SDRs are identical; that is, 

•7"PI-SDR = -^BC-SDR = -^VA-SDR- 

The proof will be described in the next two subsections. From Theorem Q] we make the important 
conclusion that 

Corollary 1 For 16-QAM MIMO detection, the relaxation problems of PI-SDR, BC-SDR, and VA-SDR 
[given in ©, (0, and ©, respectively] are equivalent. In particular, 

1) if (U,u, S,s) is an optimal solution of PI-SDR, then (S,s) is an optimal solution of BC-SDR; 

2) i/(B, b) is an optimal solution of VA-SDR, then (WBW T , Wb) is an optimal solution of BC-SDR; 

3) if (S*, s*) is an optimal solution of BC-SDR, then there exists (U*, u*) such that (U*, u*, S*, s*) 
is an optimal solution of PI-SDR; and 

4) if (S*, s*) is an optimal solution of BC-SDR, then there exists (B*, b*) such that (WB*W T , Wb*) = 
(S*,s*) and (B*,b*) is an optimal solution of VA-SDR. 

Some further discussions are as follows. 

1) The SDR equivalence established here does not depend on the objective function /. Hence, the 
equivalence applies also to other similar SDR applications where the objective function is different. 
One such application is blind ML detection of orthogonal space-time block codes [22], in which 
/ is quasi-convex and is rather different from what is considered in this paper. 

2) From Corollary Q] we see that an optimal BC-SDR solution can be directly obtained from an 
optimal PI-SDR or VA-SDR solution. In fact, an optimal PI-SDR or VA-SDR solution can also 
be constructed from an optimal BC-SDR solution in a closed-form manner. This is because our 
proof of Theorem Q] is constructive, which shows how a BC-SDR feasible point can be converted 
to a PI-SDR or VA-SDR feasible point (and vice versa). Simply speaking, constructing a PI-SDR 
solution is rather straightforward, and constructing a VA-SDR solution is more involved requiring 
some matrix factorization and vector decomposition procedures; for the construction details readers 
are referred to the proof in the following subsections. 
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3) The three SDRs can be proven to be equivalent for larger QAM sizes. For the equivalence of VA- 
SDR and BC-SDR, the proof can be generalized using a similar principle. But, for the equivalence 
of PI-SDR and BC-SDR, the proof is much harder and tedious even for 64-QAM. This will be 
elaborated upon in the next section. 

The proof of Theorem 1 consists of two parts: proving that ^bc-SDR = -^pi-SDR. and -Fbc-SDR = 

•7"VA-SDR- 

B. First Part of the Proof of Theorem 1: J^bc-SDR = -^pi-SDR 

We first show that if (U, u, S, s) G Wpi_sdr, then (S, s) is feasible to .Fbc-SDR- Given (U, u, S, s) G 
Wpi-sdr> the PI-SDR feasibility condition U y uu T implies that Uu > uf for alii = 1, . . . , N. Hence, 

= Uu - Win + 9 > Ui - lOui + 9 = (ui- \){ui - 9) 

for all i. The inequality above is the same as (Su — l)(Su — 9) < 0, or 1 < Sa < 9. This shows that 
(S,s) G J^BC-SDR- 

Next, we show that for any (S,s) G .^bc-SDR, we can explicitly construct a (U,u) such that 
(U, u, S,s) G VVpi_sdr- Consider the following construction from (S,s) G ^bc-SDR : 

u = d(S) 

U = uu T + D(w) (18) 

where D : W N — > W NxN is the operator that outputs a diagonal matrix with its main diagonals being the 
input, and w is given by 

Wi = -(Su - l)(S u - 9) = -{ Ui - l)( Ui - 9), (19) 

for i = 1, . . . , N. Since 1 < Su < 9, we have Wi > 0. It follows that U — uu T = L>(w) y 0. Moreover, 
from (fT8T)-(fT9l), one can see that 

Uu - I0m + 9 = Wi + - 10ui + 9 = 

for all i. This proves that (U, u, S,s) is feasible to Wpi_sdr- 

The proof above indicates that whatever a point is feasible to ^pi-sdr it is also feasible to .Fbc-SDR, 
and vice versa. We therefore conclude that ^bc-SDR = -^pi-SDR- 
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C. Second Part of the Proof of Theorem 1: -Ta/a-SDR = -^BC-SDR 
Let X G S> N+1 and Y G S 2N+1 be two PSD matrices taking the form 



X 



S s 



s T 1 



>= 0, 



B b 

b T 1 



>- 



where (S, s) G E> N x R N , (B, b) G 8 2N x R 2N . By Schur complement, the two matrices satisfy S y ss T 
and B >- bb T . We assume (S,s) = (WBW T ,Wb), and this condition can be expressed in a matrix 
form 



X = T YT 



(20) 



where 



W T 
1 

Since Y y 0, Y can always be represented in a square-root factorization form 

Y = R T R 

for some square root factor R = [ rj , . . . , r2N+i ] 6 ]R( 2Ar + 1 ) x ( 2Af + 1 ), with ||r2;v+i|| = 1 (owing to 
||rj|| 2 = Ya and ^2jv+i,2JV+i = !)• Similarly, X can be characterized as 

X = Z T Z 

for some square root factor Z = [ z x , . . . , zjv+i ] G R( 2Ar + 1 ) x ( Ar+1 ), ||zjv+i|| = 1. We see that d20j) holds 
if 

Z = RT. (21) 



Let us partition 



R 



U 



V 



T2AT+1 



where U, V G R( 27V + 1 ) xAr . Substituting ([22]) into (Hfl, we show that Z 
equivalently 

Zi = Ui + 2vi, i = l,...,N, 

ZAT+l = T2AT+1, 



(22) 

[ U + 2V | r 2N+1 }, or 

(23) 
(24) 



where Uj is the ith column of U, and v, is defined in a similar way. 
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Now, suppose (B,b) G Sva-SDR- Since cZj(B) = 1 for all i (where dj(-) means that di(A) = An), we 
have ||uj|| = Y ti = di(B) = 1 and ||vj|| = Y i+N)i+N = d i+N (B) = 1 for i = 1, . . . , N. With (I23t-(f24b 
satisfied, it holds true that 

||zj|| < ll u ill + 2 ll v ill = 3 > 
||zi|| > 2||vj|| - ||uj|| = 1, 

for i = 1, . . . , N. This translates into an S that satisfies di(S) = Xn = ||zj|| 2 G [1, 9]. And this further 
implies that (S,s) G ^bc-SDR- O n the other hand, suppose (S,s) G ^bc-SDR- There is no problem for 
(1241) to be satisfied, and we find (uj,Vj) satisfying d23l ) by resorting to the following lemma: 

Lemma 1 Let z G K n , n>2 be a given vector satisfying 

(3 -a< ||z|| < /? + a 
rame a, /3 > 0. T/ien there exist two unit 2-norm vectors u and v smc/i that 

z = au + /3v. 

The proof of Lemma [T] is given in Appendix lAl Essentially, the proof shows how to construct (u, v) from 
z in a closed-form manner. Applying Lemma Q] to d23l) (with a = 1 and /? = 2), for each i we obtain 
(ujjVj) that satisfies d23l ) for any ||zj|| G [1,3] (or dj(S) G [1,9]) and then achieves = ||vj|| = 1 
at the same time. This means that the resultant R [cf., Eq. (l22l )1 has unit 2-norm columns, and as a 
consequence dj(B) = Y„ = ||zj|| 2 = 1. Hence, we have (B,b) G Sva-SDR- 
We have shown by construction that Ja/a-SDR = -^BC-SDR- 

V. Generalizations to Larger QAM Sizes 

Our endeavor now turns to more challenging cases in which the QAM size is larger than 16. In what 
follows, we will prove that i) for any 4 9 -QAM constellation (where q > 1), VA-SDR is equivalent to 
BC-SDR; and that ii) for the 64-QAM constellation, PI-SDR is equivalent to BC-SDR. Details regarding 
i) and ii) will be described in the first and second subsections, respectively. 

A. Equivalence of VA-SDR and BC-SDR for A q -QAM 
For 4 9 -QAM, the ML problem to be addressed is 
min ||y — Hs|| 2 

s.t. Si G {±1,±3,±5, ... ,±(2« - 1)}, i = l,...,N. 
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Its virtually antipodal formulation takes the form 

min ||y-HWb|| 2 

be{±i}" N 

where 

W = [ I 21 41 81 . . . 2 q - x \ } G R Nx « N 



and 



with hi G M. N for all i. Again, both the VA-SDR and BC-SDR problems in this case can be represented 
by the expression 

min f(S,s) 

(S,B)£^ 

where the feasible set T for BC-SDR is defined as 

^bc-Sdr = { (S, s) | S y ss T , 1 N < d(S) ■< {21 - lfl N } (25) 

and the feasible set for VA-SDR is 

^VA-SDR = { (S, s) = (WBW T , Wb) I (B, b) G £>VA— SDR } (26) 
tfvA-SDR = { (B, b) | B y bb T , d(B) = l gN } (27) 

It is shown that the equivalence of BC-SDR and VA-SDR is promised even for higher-order QAM. 

Theorem 2 Consider a 4F-QAM constellation, where q > 1. It holds true that 

-^VA-SDR = -^BC-SDR- 

The proof of Theorem |2] is given in Appendix IB1 It is a generalization of its 16-QAM counterpart in 
Section llV-CI Like the 16-QAM case, the proof reveals the possibility that an optimal BC-SDR solution 
can be used to construct an optimal VA-SDR solution in an analytical fashion, or vice versa. 

B. Equivalence of PI-SDR and BC-SDR for 64-QAM 

The original work of PI-SDR [8] concentrates only on the 16-QAM constellation, but it is clear from 
that work that the idea can be extended to the 64-QAM constellation. To see this, we start with the 
following 64-QAM ML formulation 

min ||y — Hs|| 2 

s.t. sj G {ri,r 2 ,r 3 ,r 4 }, i=l,...,N 



SUBMITTED TO IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING, AUG 2008 



14 



where {r±, r 2 , r^, r 4 } = {1, 3 2 , 5 2 , 7 2 }. The idea is to consider the polynomial characterization 

4 5 

u G {ri,r 2 ,r 3 ,r 4 } JJ(u - n) = ^peu 1 ' 1 = 

i=l £=1 

where {pj} is the set of polynomial coefficients associated with the roots {r{\. Like the development in 
16-QAM PI-SDR, we reformulate the ML problem as 

min f(S,s) 
U,u,S,s V ; 



S.t. 



U 



Uu 


Ul 2 


,u = 


Ul 


uf 2 


u 22 _ 




u 2 



(28) 



S = ss T ,U = uu T 
d(S)=ui,d(U u )=u 2 

p x l N +p 2 ui +p 3 d(Uu) +Pid(\Ji 2 ) +P5d(V 2 2) = 0. 
The formation in (|28T ) is valid because its constraints essentially restrict m = sf, dj(Uu) = uf, di(XJi 2 ) 
uf, and d;(U 22 ) = uf. From <2g), we obtain the 64-QAM PI-SDR: 

min f(S,s) 

U,u,S,s V ^ 



S.t. 



U 



Uu 


Ul 2 




Ul 






,u = 




.Uia 


u 22 _ 




u 2 



S y ss T ,U ^ uu T 
d(S) = ui,d(Un) = u 2 

Pll N + p 2 U 1 +P3^(Un) + J» 4 d(Ui 2 ) +P5^(U22) = 0. 

For BC-SDR, the relaxation is given by 

min /(S,s) 

s.t. S y ss T , nl N ^ d(S) ^ r 4 ljv- 
The main result here is presented as follows: 



(29) 



(30) 



Theorem 3 Consider a general situation where the roots {r^} are allowed to be arbitrary (not necessarily 
the roots in 64-QAM), and assume < r\ < . . . < r 4 < oo. The PI-SDR problem in j2% and the BC-SDR 
problem in (1301) are equivalent in yielding the same feasible set corresponding to (S,s) (and thus the 
same optimal solutions), under the following sufficient and necessary condition 



y/r 4 - r~i < mm{^r 3 - r[ + \Jr 2 - ri, V r A~ r 2 + V r 4 - r^}. 



(31) 



It can be verified that the 64-QAM roots ({ri, r 2 , r 3 , r 4 } = {1, 3 2 , 5 2 , 7 2 }) satisfy (|3"T1 ). We therefore 
conclude that 
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Corollary 2 For the 64-QAM constellation, the PI-SDR problem in d29t and the BC-SDR problem in 
(1301) are equivalent in yielding the same feasible set corresponding to (S, s). 

Proof of Theorem\3} The proof is far from trivial compared to its 16-QAM counterpart. Consider the 
following lemma shown in Appendix ICl 



mm 

V 1 ,...,V JV ,S,s 



V,; 



(32) 



^0, i = l,...,N 



Lemma 2 The PI-SDR problem in d29b is equivalent to the following alternate PI-SDR problem 

/(S,s) 

s.t. S y ss T ,d(S) = [ v u , . . . ,Ujv,i } T 

«il V i2 V i3 
Vi2 V i3 

Pi + Yd=\ Vt+iVii = 0, i = l,...,N 
in the sense that the feasible sets corresponding to (S,s) are identical for the two problems. 

The proof of Lemma [2] follows the same approach as the equivalence proof for the 16-QAM PI-SDR 
and BC-SDR (in Section IIV-BI ). However, by Lemma |2] alone, we are unable to see the equivalence of 
the 64-QAM PI-SDR and BC-SDR immediately. 

To gain further insights, let us re-express the alternate PI-SDR formulation in d32l as 

min /(S,s) 

s.t. S t ss T , Sa G V, i = l,...,N 

where we define 



V = { S G E | S = [V]i2, V G V } 

V = i V G S 3 



|vg§ 3 V^0,V = Hank((l,v)),pi +^p t+1 v t = 0,v G M 4 ] 



(33) 
(34) 



with the operator Hank : 



p2n-l 



standing for 



Hank(ai, . . . ,a 2n -i) 



ai a 2 ... a n 
a 2 ' • ■ : 

: •■ a2n~2 

On ... <32n-2 ^2n-l_ 

Our interest now turns to analyzing the set V, which has to be done by analyzing V. Consider the 
following lemma proven in Appendix |Dj 
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Lemma 3 The set V in (1341 ) is equivalent to 

1 



(35) 



i=i i=i ' 

where &£ = [ 1 r$ rf ] T . 

Lemma [3] provides an interesting implication. To describe it, let 

conv{aiaf , . . . , a 4 a 4 } = i V V = ^ 9 e a. e aJ , h 0, ^ 9 e = 1 I 

^ i=\ t=\ > 

be the convex hull of {sliSlJ , . . . , a 4 aj }. It can be verified from (1331 ) that V 5 conv{a!a^, . . . , a 4 aj }, 
though V C convjaia^, . . . , a 4 aj } is generally not trudj. Consequently, we have 



V D { S = [V]i2 | V G conv{ ai a{ , . . . ,a 4 ai} } 

= ! S = J2^i4]l2 0>rO,£> = lj 
^ €=i £=i ' 

= { S = E^ »>r 0,^^ = 1 } = [n,r 4 ]. 

£=1 £=1 ' 

T/iw implies that the 64-QAM PI-SDR is no tighter than the 64-QAM BC-SDR. But we also show in 

Appendix |E] that 

Lemma 4 Lef < r\ < . . . < r 4 < oo. We have D = [r±, r 4 ] if and only if 



\A"4 - H < min{\/?"3 - ri + V r 2 - ri, V r 4 - ^2 + - ^3}- 

As a result, PI-SDR can be equivalent to BC-SDR under the condition in Lemma 0] thereby completing 
the proof of Theorem [3] 

VI. Some Further Simulation Results 

We provide two more sets of simulation results. In the first set, the purpose is to verify the SDR 
equivalence for the 64-QAM and 256-QAM cases. The simulation settings are the same as those of 
the 16-QAM simulation example in Section UlLDl and the MIMO size is (M,N) = (4,4). Simple 
rounding is employed for the SDR methods. The results are plotted in Fig. [2J We see that the symbol 
error rates (SERs) of the PI-SDR, BC-SDR, and VA-SDR with simple rounding II are generally identical, 
which corroborates our theoretical results. It is also noticed that the performance of 64-QAM PI-SDR 



3 By numerical test, it was found that there exists a V £ V such that some of the constituent Bi can be negative. 
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slightly deviates from that of 64-QAM BC-SDR and VA-SDR at SNR= 45dB. We found that this was 
due to some numerical problems encountered by the interior-point SDP solver (which is SeDuMi [27] 
here). In fact, the polynomial coefficients in 64-QAM PI-SDR have values ranging from p§ = 1 to 
pi = 1 x 3 2 x 5 2 x 7 2 = 11,025. Such a large dynamic range could be the cause of the numerical 
inaccuracy. Moreover, Fig. [2] illustrates that VA-SDR with simple rounding I is not working in 64-QAM 
and 256-QAM (see Section ITlI-D I for the definition of simple roundings I and II). This problem, which has 
also been noticed by Mao et al. [10], may partially be answered by the equivalence proof for BC-SDR 
and VA-SDR; cf., Section ITV-CI and Appendix B. In essence, the derivations there revealed that the VA- 
SDR solution with respect to (B,b) may be non-unique, even though its BC-SDR counterpart [in form 
of (S,s)] is unique. In particular, a key component, namely Lemma [T] is not a unique decomposition. 




256-QAM VA-SDR (Simple rounding I) 



256-QAM VA-SDR (Simple rounding II) 

10 

5 10 15 20 25 30 35 40 45 

SNR, in dB 

Fig. 2. Comparison of symbol error rates of PI-SDR, BC-SDR, and VA-SDR in an 4 x 4 system with either 64-QAM or 
256-QAM. 

In the second set of simulation results, a comparison is made between SDR and some other MIMO 
detectors for the 64-QAM case. Those detectors include the zero-forcing (ZF) detector, the optimal sphere 
decoder, the lattice reduction aided ZF (LRA-ZF) detector [29], [30]. It is worthwhile to mention that 
the LRA-ZF detector has been shown to achieve the full receive diversity [31]. We tested the BC-SDR 
method only, as the other two SDR methods will provide similar results. A specialized interior-point 
SDP solver [21] is used to handle the BC-SDR optimization. For solution rounding, we employ a better 
procedure, namely the Gaussian randomized rounding. The procedure has been described in [8], and the 
number of randomizations is set to 100. The results are plotted in Fig. [3] (a similar set of results for 
16-QAM is also available in [21]). Fig. [3ja) shows the case of M = TV = 8, where we see that LRA-ZF 
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->fr- LRA-ZF \ 
No-intcrfcrcncc lower bound 

-5 5 10 15 20 25 30 35 40 

SNR, in dB 
(c) (if, N) = (40, 40) 

Fig. 3. Comparison of the various detectors in 64-QAM systems. 



detector gives better performance than BC-SDR except for some rather low SNR points. We increase the 
problem size to M = N = 16 in Fig. 0b). For this problem size it is computationally too hard to run 
the optimal sphere decoder. The figure illustrates that for SNRs less than 27dB, BC-SDR outperforms 
LRA-ZF. Let us further increase the problem size to M = N = 40. As illustrated in Fig. He), now 
BC-SDR exhibits further improved performance compared to LRA-ZF. 

The comparisons above somehow suggests that SDR has significant performance advantages for large 
problem sizes and/or for low to moderate SNRs. 
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VII. Conclusion and Discussion 

This paper analyzes the relationships of three SDR-based MIMO detection methods for high-order 
QAM, namely PI-SDR, BC-SDR, and VA-SDR. We have proven that the three SDRs are actually 
equivalent, despite their different appearances. The essence of the equivalence is that an optimal solution 
of one SDR can always be constructed from that of another SDR. The equivalence result is not restricted 
to MIMO detection. It can be applied to other problems such as blind ML detection of orthogonal 
space-time block codes [22]. 

Since the three SDRs are now known to be equivalent, the comparison should turn to their computational 
costs. One can check numerically that BC-SDR is the cheapest computationally, followed by VA-SDR and 
then PI-SDR. Hence, it sounds that BC-SDR should be the method of choice among the three methods. 
While this is our present conclusion, our opinion is that each of the three SDRs is interesting in its own 
right where a different way of utilizing QAM structures has been demonstrated. Moreover, the exposition 
of the PI-SDR and VA-SDR ideas might help inspire future works for devising stronger SDR methods. 

In the introduction we have mentioned the work by Yang et al. [11], who proposed an improved 16- 
QAM BC-SDR method by tightening the feasible set with more constraints. That work showed that the 
tightened BC-SDR method can perform better than the 16-QAM PI-SDR. Using the equivalence result 
here, one can actually infer that the tightened BC-SDR can perform better than the 16-QAM VA-SDR 
as well. But, presently it is not known if there is any connection of the above-described SDR methods 
to the SDR work by Mobasher et al. [12]. Studying this would be an interesting future direction. 

Appendix 

A. Proof of Lemma [7J 

The proof is constructive. Let 

u = e-^-r + \/i - e 2 z ± (36) 

|z| 

where z± is a unit 2-norm vector orthogonal to z, and \9\ < 1. It can be verified that the u in (l36l ) 
satisfies ||u|| = 1. Moreover, let 

v = -(z-em) (37) 

which is purposely constructed to satisfy z = au + /3v. One can show from (l37l) that the unit norm 
condition llvll = 1 is achieved when we choose 
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Now the remaining problem is the condition under which \9\ < 1. It can be verified from (|38T ) that if 

(3 — a < ||z|| < (3 + a then \9\ < 1 is guaranteed. 

B. Proof of Theorem [2] 

The idea is similar to the proof in the 16-QAM case, described in Section ITV-CI We consider two PSD 
matrices 



X 



S s 



s T 1 



>z 0, 



B b 

b T 1 



>- 



that satisfy (S,b) = (WBW T , Wb). That condition is shown to be achievable if 

W T 



Z = R 



1 



(39) 



where Z G R(</W+i)x(iV+i) and R £ R (qN+i)x( q N+i) are square root factors of x and Y, respectively 
(or Z T Z = X, R T R = Y). The objective is to show that if (S,s) G .^bc-SDR* then we can construct a 
(B,b) G £>va-SDR satisfying ( f39b ; and vice versa. 
For general 4^-QAM where W is expanded to [ I 21 . . . 2 9_1 I ], Eq. d39l can be rewritten as 

9-1 

Zi = J2 2jr i+jN, i = l,...,N, (40) 



3=0 

%N+l = XqN+1- 



(41) 



To achieve (|4TT ) is easy, and the nontrivial part lies in (l40l . First, suppose (B,b) G £>va-SDR- Then the 
resultant R satisfies ||r-j|| = 1. The vectors Zj satisfying (l40l would then have bounds 



g-l g-1 

N| < ^^'Uri+jjvll = Y, 2 ' = 29 ~ h 

3=0 3=0 



|zi|| > 2? ll r i+(g-l)Ar| 



9-2 



9-2 

^2 J r i+j7V 

j=0 



> 2 9 " 1 - 2 J ' = 2 9 " 1 - (2 9 ' 1 - 1) = 1. 

j=Q 

This means that the corresponding S has dj(S) = ||zj|| 2 G [1, (2 9 — l) 2 ]. Hence, (S,s) G ^bc-SDR- 
Second, suppose (S,s) G -Sbc-sdr- Let us choose, for each i = 1, . . . , N, 



ri = r i+N 



'i+(q-2)N 
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for some Uj,Vj G IR' 7iV+1 . The condition in (l40l) becomes 

q-2 



(42) 



Using Lemma [U we can construct a (iij, Vj) that satisfies (|42]) for any ||z;|| = -y/dj(S) G [1, 2 q — 1] while 
achieving ||uj|| = ||vj|| = 1. The resultant B will therefore satisfy cij(B) = Ya = 1 for all z, meaning 
that (B,b) G Sva-SDR- 

The proof of Theorem [2] is complete. 



C. Proof of Lemma\2\ 

Suppose that (U, u, S, s) is feasible to the original 64-QAM PI-SDR problem in ([29]). Set 



V,: 



1 



e 





t 
e i 

T 

i+N 



1 


T 

U 




1 








u 


U 







e, 





(43) 



for i = 1, . . . , N, where G M? N is a unit vector with the nonzero element at the ith element. It follows 
from U y uu T and (03]) that Vj y for all i. Moreover, (03]) equals 

1 Uu U 2 i 

Vi= u u [Un]« [Ui2]« • (44) 

U2i [Ul 2 ]ii [U 22 ]ii 

Since = [Un]ii, every Vj in (1441 satisfies the Hankel structure in the alternate PI-SDR problem in 
(l32l) . It also follows from (l44b that the equality constraints arising from the polynomials are satisfied. 
Hence, (Vi, . . . , Vjv, S, s) is feasible to the alternate PI-SDR problem in d32l ). 

On the other hand, suppose that (Vi, . . . , V^, S, s) is feasible to the alternate PI-SDR. Set 

ux = [v 11 ,...,v N>1 ] T , (45a) 

Uu = C 2 - D(ui ui) + uiuf , (45b) 

u 2 = d(Uu) = [ vu, ■ ■ ■ ,vn,2 ] T , (45c) 

Ui 2 = C 3 - D(ui u 2 ) + uiuf, (45d) 

U 22 = C 4 - L>(u 2 u 2 ) + u 2 u^, (45e) 

where is the Hadamard product, and 

Cj = D(yij, . . -,v N)i ). 
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It can be shown that (1451) satisfies the equality constraints of the polynomials. Let us examine if 
U >- uu T . We see that 



U- uu 



T 



C a -D(ui©ui) C 3 -D(ui0u 2 ) 
C 3 - D{\ii u 2 ) C 4 - D(u 2 u 2 ) 
It can be shown using basic matrix results that the specially structured matrix above is PSD if and only 
if 

,2 



(46) 



Vi,3 ~ Ul,iU2,i V it 4 - 

are PSD for all i = 1, . . . , N. By using = v^i and U2,i = t>i. 2 and Schur complement, it is shown that 
(l46l ) are indeed PSD. We therefore conclude that (Vi, . . . , Vjv, S, s) is feasible to the alternate PI-SDR 
problem. 



D. Proof of Lemma\3\ 

Consider a problem setting as follows: Let {ri, . . . ,rif\ be a given set of roots, and assume that the 
roots are distinct. Consider the following two sets 

Vi = | V = Hank(v) v G R L+1 , v x = 1, p T v = 0, V h 

where p = [ pi, . . . ,pl+i } T contains the polynomial coefficients corresponding to {ri, . . . ,r^}, i.e., 
J2e=i Per 1 " 1 = for all r G {ri, . . . , r^}; and 



V 2 



e=i t=i ' 



where = [ 1 r£ r| . . . r^ 2 ] G M L / 2+1 , and L is assumed to be even. Our objective is to prove that 
Vi = V 2 - Clearly, Lemma [3] is a special case of the above problem where L = 4. 

By definition, every V G Vi can be parameterized by some v G such that v\ = 1 and p T v = 0. 

Let 

b t = [ l,r e ,r 2 e ,...,r£ ]ei L+1 

for t = 1, . . . ,L. Since every contains one of the true roots, it satisfies p T b^ = 0. Hence, we have 
the following condition to satisfy 



p T [b!...b L v]=0. 



(47) 
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The submatrix [ bi . . . t>£ ] G r( l + 1 ) xL is linearly independent, being Vandemonde with distinct roots. 
Subsequently, (|47T ) can be satisfied only when 

L 

1=1 

for some coefficients 6 G R L . Since 1 = v\ = ^2d = \6^i]i = J2e=i ®h tne coefficients satisfy 
X^=i &i = 1- Moreover, by noticing that 

Hank(b £ ) = a e aj , (48) 

we have 

L 

V = Hank(v) = ^ 9 e a e aJ . 

Hence, any V G Vi lies in V2. 

Likewise, it can be verified that every V G V2 lies in V\. For every V G V2 which can be characterized 

as V = J2e=i ^a^aj, J2e=i Ql = 1» set 

L 

1=1 

It follows from (|48T ) that Hank(v) = V. Moreover, this v satisfies v\ = J2e=i @t\bt]l = 1» anc ^ P Tv = 
Eli OeP T b e = 0. 

£. Proof of Lemma |4] 

By Lemma [3l the set I? can be expressed as 

1 



-44 4 v 



This set is a closed convex set, and therefore must be in form of an interval [L, U}. The proof is divided 
into three parts: solving L, solving U, and integrating the results. 

1 ) Solving the lower bound: We find the lower bound L by solving the problem 

4 

L = min 

" 1=1 

4 4 

s.t. ^ O^aJ ^ 0, = 1 





£=l (49) 



=1 1=1 

^3 



Let Xj = 0j.fi, i = 1, 2, 3. Using 0i = 1 — ^ i=1 £j, Problem d49l can be re-expressed as 

3 

L = min n + y^Xj(rj + i — 7*1) 

X * ' 

1=1 3 (50) 
s.t. aiaf + ^2 Xi(a i+1 aJ +1 - a x af ) y 0. 



i=l 
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By strong duality, solving (1501) is the same as solving its dual which can be shown to be 



L = max n — tr(aiaTZ) 

s.t. z y o, (51) 

tr[(aiaf - a i+1 af +1 )Z} = r\ - r i+x , i = 1,2, 3. 
From the objective in (|5T| ), it is clear that L = n if and only if Z is feasible and satisfies 

tr(aiaf Z) = 0. (52) 

Let us consider the construction of a PSD Z satisfying (l52l . Eq. (l52l implies that Z has rank no greater 
than 2. Thus any such PSD Z can be represented by 

Z = RR T 



where R G Ir is such that R a! = 0. Such an R can be parameterized as 



R= [ Woi, Wo 2 ], 



for some 0:1,0:2 G M 2 , where 



W 



1 1 

■l/n -2/n 

l/rf 



(One can easily check that af W = 0, thereby R T ai = 0.) Therefore, any PSD Z satisfying (l52l can 
be expressed as 

Z = WGW T 



(53) 



where 



G 



a b 
b c 



01 02 



a; 



y 



can be any 2x2 PSD matrix. 

By substituting the matrix form in d53l into the equality constraints in d5T| ), we obtain 

af +1 WGW r a !+1 = r l+1 - n, i = 1,2,3. 



(54) 



We seek to find the sufficient and necessary conditions for satisfying (I54I ). By noticing that 

1 

1 _i I*.. n 

W r a m 



1 -i/n 

1 -2/n l/rf 



1 - r i+ x/r\ 
(1 - r m /n) 2 
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Eq. (l54l ) can be decomposed to 

ri - a(r i+1 /n - 1) + 26(r m /n - l) 2 - c(r i+1 / n - l) 3 = 0, (55) 
for i = 1,2,3. Let us define a polynomial function 

f(u) = r\ — au + 2bu 2 — cm 3 . 
Since the function satisfies /(rj+i/ri — 1) = for i = 1, 2, 3 [cf., Eq. (|55T)1. it permits a factored form 



3 

f(u) = n - n/u j; , 

where 



'y+l/ri - 1 > 0. 



= r 

By expanding the factored form of f{u) to the polynomial form, we determine (rather tediously) that 

a = {u\U2 + U2U3 + > 0, 

U1M2W3 

-(Ui + U 2 + U 3 ) > 0, C= > 0. 



2uiU 2 U3 ' U1U2U3 

The remaining part lies in ensuring that the resultant G is PSD. We already have a > and c > 0, so 
the last condition is b 2 — ac < by Schur complement. With some cumbersome derivations, we show 
that 

b 2 -ac= ( — — ) [U 3 - {y/u x - y/u 2 ) 2 ][uz - (y/u 1 + Vu 2 ) 2 }. 

In order to achieve b 2 — ac < 0, we need 

- ^2? <u 3 < (Vui + V^2?- (56) 



Summarizing, we have L = r\ if and only if (1561 ) holds. 

2) Solving the upper bound: The method of the proof is exactly the same as the previous, and hence 
the detailed derivations are omitted for brevity. Essentially, we consider solving the upper bound 

U = max Yh-i ^i r i 
e 

s.t. Eti^af^O, Eti^ = l 

by solving its dual 

U = min ta + trfa^aTZ) 

zes 3 

s.t. z y 0, (57) 
tr[(a 4 aj - a^af )Z] = r { - r 4 , i = 1, 2, 3. 
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From (I57T ) it is shown that U = r 4 if and only if 

{V^2-V^3? <Vi < (^ 2 + ^ 3 ) 2 - (58) 

where vi = 1 — rj/r 4 > for i = 1, 2, 3. 

3) Combining the conditions: The final task is to combine the conditions in d56l ) and d58l ). We can 
express (l56l) as 

\A"3 - - V r 2 - r% < V r 4 - ri < V r 3~ ri + V r 2 - H- 
The lower bound is redundant because for any r 4 > . . . > n > 0, 

V^4 - n > V^3 - n > \/r 3 - ri - \/r 2 - r x . 
Moreover, d58l) can be expressed as 

y/r^ - r 2 - \/r 4 - r 3 < \/r 4 - ri < \/r 4 — r 2 + \/r 4 - r 3 , 

and again the lower bound can be shown to be automatically satisfied. We therefore obtain the sufficient 
and necessary condition in Lemma 01 completing the proof. 
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